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The recent identification of a time and temperature dependent spin glass correlation length, £(t w ,T), 
has consequences for samples of finite size. Qualitative arguments are given on this basis for depar- 
tures from t/t m scaling for the time decay of the thermoremanent magnetization, MrRM(t,t w ,T), 
where t is the measurement time after a "waiting time" t w , and for the imaginary part of the ac 
susceptibility, Consistency is obtained for a more rapid decay of Mtrm {t, t w , T) with 

increasing t w when plotted as a function of t/t w , for the deviation of the characteristic time for 
MrRM{t,t w ,T) from a linear dependence upon H 2 at larger values of H, and for the deviation of 
the decay of x"(w,t) from cut scaling upon a change in magnetic field at large values of wt. These 
departures from scaling can, in principle, be used to extract the particle size distribution for a given 
spin glass sample. 

PACS numbers: 75.50.Lk, 75.10.Nr, 75.40. Gb 
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I. INTRODUCTION 

Slow spin-glass dynamics have been investigated 
through examination of the time dependence of the ir- 
reversible magnetizationJ!l aad the out-of-phase dynam- 
ical magnetic susceptibility.cl Both are related, but typt 
ically explore different time domains. A cecent papertl 
has shown how the Parisi order parameter ,□ x(q), can be 
extracted from high precision measurements of the decay 
of the thermoremanent magnetization, MTRAi(t,t w ,T), 
where t is the measurement time after waiting a time 
t w at the measurement temperature T(< T g ). The time 
scale of the experiments restricts the examination to a 
small region of overlap space, but the shape and magni- 
tude of the results are in quantitative accord with mean 
field expressions. □ 

The barrier model of NemotcO and Vertechi and 
VirasordU has been adapted to the metastable states ob- 
served in experiment Dynamics are ascribed to bar- 
rier hopping within ultrametric manifolds of constant 
magnetization Ma Measurements by Chu et aLxSl along 
with other observations of a similar nature,EU suggest 
that a change in magnetic field reduces the barriers in 
the initially occupied magnetization manifold. This re- 
duction can be conceptually thought of as diffusion be- 
tween states of constant J\d through intermediate states 
of lower Zeeman energy.EJ We designate the reduction 
in each of the barrier heights by a Zeeman energy, E Zl 
proportional to the number of participating spins, N s , 
which lie within a coherence length £ (t w , T) from one 
another. This length scale has been extracted from mea- 
surements on both insulating and metallic spin glasses.E3 
As shown in Ref. 12, and exhibited below in this paper 



as a consequence of further measurements, the spin glass 
correlation length appears universal, having the value 
£(t w ,T) = (t w /T ) aT/T *, with a = 0.153, for all spin 
glasses measured so far (Cu : Mn, Ag : Mn, and the 
thiospinel CdCri/jIno.sS/i). The length £(t w ,T) is in 
units of the typical spin-spin spatial separation. 

This paper examines the consequences of a finite length 
scale associated with spin glass order. In particular, we 
address the question of what happens when £(t w , T) ex- 
ceeds a physical length r associated with finite sample 
size {e.g. defects in an insulating structure, or crystal- 
lites in a polycrystalline sample). Typical laboratory 
waiting times and measurement temperatures result in 
!;(t w , T) ss 10 — 100 nm, approximating the coarseness of 
powdered or polycrystalline samples. Such materials do 
not possess a single grain size r, but rather a distribution 
of grains sizes, P(r). This paper will show that it may be 
possible to extract P(r) from time dependent magnetic 
measurements . 

A consequence of £ (t w , T) w r is that conventional 
scaling relationships may be violated. This paper will 
associate the lack of scaling with t/t w found-, Jar the 
time decay of MrRM{t,t w ,T) by Alba et a/JiM-il with 
£(t w ,T) ks r. Other observations which wc believe can 
be associated with £ (t w , T) m r-Ejpaj e the deviation of E z 
from a proportionality to 7? 2 J13ll3 and the deviation of 
x"(u>,t) frorp_j^i scaling upon a change in magnetic field 
for large Lut.13 

This paper will attempt to make plausible the relation- 
ship between finite size and, lack of scaling in the same 
spirit as in Bouchaud et aLt£l Sec. II describes the model 
which underlies our analysis. Sec. Ill outlines the experi- 
mental observations which require a departure from scal- 
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ing. The relationships between finite size effects and the 
lack of scaling are developed in Sec. IV. Key experiments 
and numerical simulations required to fully justify these 
relationships are discussed in Sec. V. Sec. VI presents 
our conclusions. 



II. DESCRIPTION OF THE MODEL 

The model we invoke to quantify our analysis is 
based on an extension of the "pure states" ultramet- 
ric geometry foundO for the mean field solutionQ of the 
Sherringtoni-Kirkpatrick infinite range model.E3 Previous 
experimentsS have shown that the infinite barriers which 
separate pure states derive from finite barriers separating 
metastable states, diverging at a characteristic tempera- 
ture T* as the temperature is lowered. Withm the time 
and temperature scale probed experimentally^ the value 
of T* depends only on the height of the barrier at a given 
temperature, and a universal form for dA/dT vs T was 
derived from experiment. 

Dynamics are extracted from the assumption that the 
metastable states possess the same ultrametric geome- 
try as pure states, the barriers separating the metastable 
states increasinaJinearly as the Hamming Distance D 
between statespQ A(D) oc D, and activated dynamics. 

The time decay of the thermoremanent magnetization 
follows from assuming that a temperature quench in con- 
stant magnetic field isolates the spin glass states into spe- 
cific points within phase space. This is represented by a 
probability density P(D), with D the Hamming distance 
defined by D = ^(qEA — gl_where qEA is the Edwards- 
Anderson order parameter J13 and q the overlap between 
the states separated by D. Immediately after quench, 
P(D) = 6(D). Keeping the magnetic field constant, the 
system diffuses from D = to states with D ^ as a 
function of the "waiting time" , t w , according to activated 
dynamics. As stated earlier, the ultrametric tree upon 
which this diffusion develops is one of constant magneti- 
zation Mf c associated with the field cooled stateB The 
amplitude of the delta function diminishes with increas- 
ing t w , with the associated occupancy of states at finite 
D increasing with t w . Activated dynamics in phase space 
leads to a maximum barrier surmounted in the time t w of 
magnitude A(t w ,T) — kBT£n(t~,lxp). Detailed balance, 
with experimental confirmationJjll3 leads to equilibrium 
occupation of the metastable states. .— . 

Experiments and analysis by Vincent et aZ.til establish 
that, after waiting a time t Wl cutting the magnetic field 
to zero diminishes each barrier height by an amount to 
which we shall refer as £j z , the Zeeman energy. This 
reduction was interpreted!] as the origin of the so-called 
"reversible" change in the magnetization. Their model 
assumes that all of the states occupied for D < De z [in- 
cluding those at -P(O)] immediately empty to the zero 
magnetization manifold M = 0. If A(D) is known, then 
Ds r is the value of D at which E z — A(D). Field cycling 



experiments showB that D is "respected," that is, the 
population of the states within < D < De z in the M 
manifold rapidly transitions to the states < D < De z in 
the M — manifold. At the measurement time t, where 
the "clock" starts when the magnetic field is cut to zero, 
i.e. at t w , the population of the states remaining behind 
for D > De z in the M manifold have a total magnetiza- 
tion MTRM(t, t w ,T). They decay to the M = manifold 
by diffusing from D > Ds r to the "sink" created by 
the magnetic field change D < De,- The characteristic 
response time is set by the time it takes for the popu- 
lation in the states which have surmounted the highest 
barrier, -Da(*„,t)! to diffuse to the Hamming distance at 
the edge of the sink, De z - By supposition, this is equiv- 
alent to surmounting a barrier of characteristic height 
A{t w ,T)-E z . 

For very small E z , this diffusion process yieldsEJ 
a peak in the spin glass relaxation rate S(t) — 
d[—M TRM (t,t w )/H]/d£n t at t « t w . For finite E z , 
the peak in S(t) is shifted to shorter times tffi , and was 
first noted for experiments upon tie insulating thiospinel 
CdCr1.7Ino.3S4 by Vincent et al^A and the amorphous 
system (Fe x Nii_ x )7 5 Pi 6 B 6 Al3 by Djurberg et aL.Eil The 
effective characteristic time immediately follows from ac- 
tivated dynamics, 

A{t w ,T)-E z = k B T{tntf -in t ) . (1) 

The final piece of the puzzle is the magnitude of E z . 
The overall magnetization before the magnetic field is cut 
to zero, Mf ci is essentially constant at the measurement 
temperature during the waiting time t w . One can think 
of Mf c as arising from the population of the states in the 
M manifold, each of which has the same magnetization. 
If we define the susceptibility per spin in the field cooled 
state as x/c, then the magnetization per occupant of each 
state is just XfcH . The Zeeman energy, E z , is the amount 
by which pthe barriers in the M manifold are reduced. 
BouchaudiJ ascribes this "reduction" to diffusion out of 
the constant M plane, caused by a "tilting" of the overall 
energy surface. 

For the barriers A(D) to be uniformly reduced, there 
must be a coherence associated with the "hopping" pro- 
cess. That is, there must be a certain number of spins, 
N s , which are rigidly locked together, and which partic- 
ipate in the hopping process as a coherent whole. But 
these spins possess a net magnetization N s Xf c H, or a 
Zeeman energy E z = N s Xf c H 2 - 

Hence, use of Eq. (1) in experiments as a function of 
H generates absolute values for E z , and thence for N s , 
because Xfa is known. It was this analysis which enabled 
Joh et al.jB using N s oc £(t,T) 3 , to find-£(i UI , T), and to 
compare with Monte Carlo calculationsE2l by extrapolat- 
ing to time scales ten orders of magnitude shorter than 
laboratory times. 

Measurements of logiot^f vs H 2 , from Eq. (1), should 
yield a straight line, the slope of which can be used to 
obtain an absolute value for N s (t, T), and thence £(t, T). 
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As will be shown in the next Section, the H 2 dependence 
of logxotffi is found only at the small end of the magnetic 
field range. At larger magnetic fields, the dependence on 
H veers away from quadratic to more like linear. Further, 

dynamics derived from this model should scale as — , 

173 

whereas the data quoted— by Joh et al.ill Bouchaud et 
aZ.Jlj and Vincent et aL,EJ do not. Sec. Ill (following) 
displays the results of experiments which exhibit these 
deviations, forming the basis for our subsequent analysis 
in Sec. IV where we account for these deviations on the 
basis of finite size effects, along the same lines previously 
proposed by Bouchaud et alx3 



Mn2.6 at.%, measured at SACLAY), allows one to plot 
N s vs (T/Tg)£n(t w /T ) for three different physical sys- 
tems, over a wide range of reduced temperatures, T r , 
and waiting times t w . The results are exhibited in Fig. 
2. The solid line drawn through the points, setting 
£(t w ,T) = N s z , is the relationship quoted in the Intro- 
duction, 

^,T) = (Wto)°- 153T/T * , (2) 

where the unit of length is the typical spin-spin spatial 
separation. 



III. EXPERIMENTAL EVIDENCE FOR LACK OF 
SCALING 

The previous Section showed that Eq. (1) could be 
used to obtain a quantitative value for the Zeeman energy 
E z = N B Xf c H 2 . This requires E z to scale as H 2 . Fig. 1 
reproduces Fig. 2 of Ref. 12. 
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FIG. 1. A plot of logiot^ [equivalently E z from Eq. (1)] vs 
H 2 for Cu : Mn & at.% (T/T g = 0.83, t w = 480 sec) and the 
thiospinel, CdCrT_.7Ino.3S4, (T/T g = 0.72, t w = 3410 sec) at 
fixed t w and T from Ref. 12. The dependence is linear in H 2 
for magnetic fields less than 250 and 45 G, respectively, then 
"breaks away" to a slower dependence. The scale for H 2 in 
the two plots has been adjusted to the distance from the de 
Almeida-Thouless line at the respective temperatures. 

This figure exhibits a quadrative dependence for 
logwtffl on H in the very small magnetic field change 
limit, "breaking away" to a slower dependence at a 
slightly larger field, Hi_ re ak- We shall associate Hf> r eafe 
with the smallest crystallite size. Alternatively, a lin- 
ear dependence on H can be fitted to the data on the 
thiospinefJtj at over a range from small to moderate 
H , with a deviation at very small H . 

Taking the small H region slope from Fig. 1, to- 
gether with additional data recently obtained by the au- 
thors (thiospinel, measured at UC Riverside; and Ag : 
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A=T/Tg ln(t w /x ) 

FIG. 2. A plot of N 3 on a log scale vs (T/T g )ln(t w /T ) 
for Cu : Mn 6 at.% (solid circles); thiospinel measured at 
SACLAY (open triangles); thiospinel measured at UC River- 
side (solid inverted triangles); and Ag : Mn2.6 at.% (open 
circles) . 

In a similar fashion, the decay of the thermoremanent 
magnetization MTRM(t,t w ,T) does not obey simple — 

_ tw 

scaling. Fig. 3 reproduces Fig. lb of Ref. 23, a plot of 
M TRM (t, t w ,T) t 

vs — , but with an expanded scale, lne 

M fc t w 

failure to scale with t/t w is seen clearly. 

The lack of scaling arises from two sources. The first is 
emptying of the delta function at D = with increasing 
t w , appropriate to the barrier model; or, equivalently, 
from the presence, of stationary dynamics in the model 
of Vincent et aZ.E3. The second contribution arises from 
transitions between barriers, or, equivalently from the 
non-stationary dynamics. 

In order to display the part of the decay of 

trm( , w, — ) assoc j a ^ ec [ w itfi the dynamics of barrier 

hopping, or, concomitantly, the non-stationary part of 
the magnetization decay, the data of Fig. 3 are replotted 
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in Fig. 4 with- the estimated stationary contribution of 
Vincent et alE2 subtracted from the full measured value. 

The behavior exhibited in Figs. 3 and 4 was-^scribed 
to an "ergodic" time, t erg by Bouchaud et al$3l3, associ- 
ated with the occupation of the deepest trap in a "grain" 
with finite numbers of states. We shall argue that a finite 
t erg is also responsible for the behavior exhibited in Fig. 
1. 



the magnitude of the jump in when plotted as 

a function of tot, for cot large (of the order of 1,000 or 
larger), in violation of scaling. However, for small cut, 
scaling is obeyed. These features will be shown to be 
consistent with a distribution of crystallites of finite size 
in the next Section. 
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log w (t/t w ) for Ag : Mn 2.6 at.% for t w = 300, 1000, 3000, 
10000, and 30000 seconds from Refs. 14 and 23. The lack of 
scaling with t/t w is evident at large t/t m . 
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FIG. 4. Same as Fig. 3, but with the estimated stationary 
contribution subtracted from the full measured value. The 
lack of scaling with t/t w continues to be evident at large t/t w , 

and it is clearer that — T R M ( > — ) £ Qr | on p- er wa itins times 

M fc 

decays faster than for shorter waiting times. 

In addition to the lack of t/t w scaling for the 
MTRM(t,t w ,T), there is also a departure from. scaling 
for the time dependent magnetic susceptibility.liil 

Reproducing Fig. 2 of Ref. 11 in Fig. 5 for the 
thiospinel, it is seen that there is an to dependence for 
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FIG. 5. Comparison of the effect on % °f a dc field vari- 
ation in two experiments at frequencies w = 0.1 and 1 Hz 
from Vincent et al. in Ref. 11. The field variation is applied 
at a time ti (= 350 and 35 min, respectively) such that the 
product uo ■ t\ is kept constant. The curves are plotted as a 
function of uj ■ t, and have been vertically shifted in order to 
superpose both relaxations before the field variation. At con- 
stant uj ■ t, the effect of the perturbation is seen to be stronger 
for the lowest frequency (longest t\). 

In summary, three phenomena show departures from 
the predictions of the model outlined in Sec. II: (a) the 
"break away" from the H 2 dependence of logiot^; (b) 
the lack of scaling with t/t w of the time dependence of 
the thermoremanent magnetization MxRM{t, t w ,T); and 
(c) the lack of scaling with ut of x"(cj,t) at large ust. 
These separate, but related, observations are consistent 
with a particle size distribution P(r) in the sample, such 
that the spin glass correlation length £(t w ,T) becomes 
comparable with the size of a component particle r. 



IV. FINITE SIZE EFFECTS AND THE LACK OF 
SCALING 

We examine in this Section the dynamics of a small 
spin glass particle of radius r, at a waiting time and tem- 
perature where £,(t w , T) is comparable to r. Aging ceases 
in thp-psrticle when this occurs, first noted by Bouchaud 
et alMli 

The barrier model of Sec. II has the following conse- 
quences. For increasing t w , occupied states are separated 
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by barriers which increase in height according to 
A(t w ,T) = 6Mk B T g £n £(t w ,T) . 



(3) 



Immediately after the time t w when the magnetic field 
is cut to zero, Eq. (3) specifies the maximum barrier 
height surmounted by the system. Should £(t w ,T) ~ r, 
there would be no more barriers to surmount! The oc- 
cupation of all states would be at equilibrium, and aging 
in the sense of the barrier model ceases. Of course, any 
physical system will have a distribution of particle sizes, 
P(r), so that this cessation will be "smeared" in time. 
The purpose of this Section is to explore the impact of 
£(t w ,T) « r for each of the three experimental depar- 
tures from scaling presented in Sec. III. 

(a) "Break away" from H 2 dependence of logxot^/ 

As introduced in Sec. II, a susceptibility, Xfc, can 
be associated with each spin in the field cooled state, 
resulting in E z = N s Xf c H 2 . In this way, the slope of 
the plot of E z vs H 2 generates N s , the number of spins 
locked together in a coherent state. As shown in Sec. Ill, 
the actual data do appear to scale as H 2 for very small 
magnetic field changes, but this scaling breaks down at 
slightly larger fields, Hbreak ~ 170 G for Cu : Mn 6at.% 
and Hbreak — 45 G for the thiospinel, CdCr1.7Ino.3S4, 
with each Hbreak equivalent to about 152o of the respec- 
tive de Almeida - Thouless critical field EJ An alternative 
linear dependence of E z on H does describe the data 
over the field range beginning with i^—.L—and extend- 
ing to the largest magnetic field change,E3£3 but fails at 
very small field change. At the time of the publication 
of Rcf. 12, we wrote that "We do not have a satisfactory 
explanation for this change in slope." 

We believe that the departure of the plot of E z from 
proportionality to H 2 can be understood within the bar- 
rier model of Sec. II, modified to include finite size ef- 
fects. 

Consider a spin glass particle of radius r. The number 
of correlated spins would be proportional to r 3 should 
r < £,{tw,T), but proportional to £ 3 (t w ,T) should r > 
£(t w , T). Thus, E z would be less for the smaller particles 
(r < £(t w ,T)) than for the larger particles (r > £,(t w ,T)) 
at the same value of H 2 . Further, the largest barrier over- 
come on a time scale t w would be A(r) = G.OAksTgin r 
for r < £(t w ,T) and A(t w ,T) = 6.04k B T g £n£(t w ,T) for 
r > £(t w ,T). The effective waiting time, tffi ', depends 
upon the difference A — E x from Eq. (1). For small 
magnetic field changes, this means that t^/ is larger for 
the larger particles and smaller for the smaller particles 
[£n£(t w ,T) > £nr}. 

As the magnetic field change increases, E z increases 
more rapidly for the infinite [meaning ^(t w ,T) < r] size 
particles than for the smaller [meaning £(t w ,T) > r] 
size particles. This means that the peak in S(t) = 
-M TRM (t, t W) Th 



rapidly with H for the infinite particles than for the 
smaller particles. At some value of magnetic field change 
H, t%' for the infinite and smaller particles will become 
equal. This yields an increase in apparent width for S(t), 
as observed in many experiments. For yet larger H, the 
weight of all the smaller particles dominates, slowing the 
shift of the peak of S(t) with increasing H, leading to a 
less rapid decrease of Log\§t & J' with increasing H . 

We believe this to be the origin of the "break" in the 
slope of Log w t^f versus H 2 , exhibited in Fig. 1 of Sec. 
Ill, and therefore an effect of finite size. A quantitative 
fit will require knowledge of the particle size distribu- 
tion P(r). For now, this finite size effect can explain the 
behavior of the magnetic-field dependence of the charac- 
teristic response timeOHLS within the barrier model of 
Sec. II. 



(b) Lack of — scaling for MrRM(t,t w ,T) 

tw 

The barrier model of Sec. II to describe spin glass 
dynamicsi predicts scaling as a function of — . Measure- 

tiu 



H 



din t shifts to shorter times 



ments of Ocio et a/.oEj, exhibited in Figs. 3 and 4 of 
Sec. Ill show that this is not the case in the long time 
domain t > t w . That is, Fig. 4 corrects Fig. 3 for the 

stationary contribution.0 Departure from scaling as — 

tw 

is only truly present for t > t w , as can be seen in Fig. 
4. As is clearly seen in Fig. 4, the barrier hopping or 
the non-stationary dynamics results in the relaxation of 

older systems being "faster" when plotted versus — (al- 

though, of course, when plotted versus t, the olde&nthe 
system, the slower the relaxation). Bouchaud_ef alS3 re- 
called that all of the data of Alba et aZ.EJo could be 
rescaled with the response times in the spin glass scaling 
as (t w +t) fi , with fi < 1. They noted that, from a more 
fundamental point of view, the scaling variable should 

be written as — 7-, — r~r, with r* a characteristic time 

scale. Their manuscript ascribed a physical meaning to 
t* , relating it to the finite number of available metastable 
states in real samples made of 'grains' of finite size. In 
their language, "a finite size system will eventually find 
the 'deepest trap' in its phase space, which corresponds 
to the equilibrium state.M This will take a long, but fi- 
nite time t erg ', when t w exceeds this 'ergodic' time, t erg , 
aging is 'interrupted' because the phase space has been 
faithfully probed. Beyond this time scale, conventional 
stationary dynamics resume." 

The precise origin of the ergodic time scale was "...not 
easy to discuss since we do not know-precisely what these 
'subsystems' are." Bouchaud et al$3 suggested it could 
be magnetically disconnected regions (such as grains), 
the size of which was determined by sample preparation, 
and thus temperature independent; or it could be that 
the phase space of 3c? spin-glasses is broken into mutu- 
ally inaccessible regions ("true" ergodicity breaking). We 
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argue below that finite size effects can indeed account 
for this behavitp-i-Our approach is the same as that of 
Bouchaud et aZ.IIInIS, but our arguments will be based on 
the barrier model of Sec. II. 

The reasoning follows from the time dependence of the 
spin glass correlation length. From the data in Fig. 2, 

± 0.153T/T g 

and Eq. (2), £(t w ,T) — ( — ) , where r is of the 

To 

order of an exchange time (« 10 -12 sec). This means 
that, for a given waiting time t w , the correlation length 
£(t w ,T) can be larger than a particle with size r. All 
the available metastable states in that particle would be 
occupied in thermal equilibrium, and no further aging 
would take place. The largest barrier in that particle 
has magnitude^ 2 ] A(r) = G.O^tfcsTg^nr, less than the 
largest barrier in the particles for which £(t w ,T) < r. 
For these particles [large on the length scale of £ (t w , T)] , 

A(t Wl T) = kBT£n(—j. The characteristic time for 

decay of Mt rm (t,t w , T) is proportional to A - B z S0 
Therefore, for small magnetic field changes, the charac- 
teristic decay time of the small particles, proportional to 
A(r) — r 3 \fcH 2 , is less than the decay time of the larger 
particles, proportional to A(t w ,T)—^ 3 (t W7 T)xf c H 2 . Av- 
eraged over all particles, small and large, the characteris- 
tic decay time for Mtrm^, tw, T) will be less than t w , the 
magnitude of the difference depending upon what frac- 
tion of the sample contains particles of size r < £(t w ,T), 
i.e. the particle size distribution. As t w increases, more 
of the particle sizes r will be less than £ (t w , T) , thereby 
shortening the characteristic time for Mtrm^, t w ,T) de- 
cay. The characteristic time will shorten as t w increases, 
leading to a faster decay of M(t,t w ,T) with increasing 

t w when plotted as a function of — . This is exactly 

the effect posited by Bouchaud et aZ@. In addition, the 
Zeeman energy E z for the larger particles is proportional 
to £ 3 (t w ,T) — N s . Increasing t w will increase £(t w ,T), 
causing E z to increase with increasing t w . This further 
diminishes A(t w ,T) - £ 3 {t w ,T)xf c H 2 for the lar ger par- 
ticles, adding to the reduction of the characteristic time 
for decay of Mtrm^, t w , T) with increasing t w . This ad- 
ditional contribution has also been noted by Bouchaud 
et alJB 

(c) Lack of scaling for x"(u;,i) at large tot 

The jump in x"(w,t) with the application of a mag- 
netic field is an important test for any dynamical model. 
Fig. 5 displays the frequency dependence of the change 
in when plotted against to ■ t. It is seen that 

the jump is larger, the smaller u>. The origin of this 
effect within the trap modelEJ was associated with an in- 
crease of coupling to the magnetic field, the deeper the 
trap. Within the barrier model^j this effect was associ- 
ated with an increase of coupling to the magnetic field, 
the higher the barrier. The relationship of the jump in 



x"(w,t) to these non-uniform magnetic field couplings 
arises naturally from the time dependence of the spin- 
glass correlation length [Eq. (2)], £(£, T). The smaller 
ui, the larger t, when ^"(o;,t) is plotted as a function 
of oj ■ t. But larger t means larger £(£, T), thence larger 
N s [ oc £ 3 (i,T)], thence larger E z ( oc N s Xf c H 2 ). This 
increase in E z with increasing t maps directly onto thn 
non-uniform magnetic £f Id coupling of the trap modelEll 
and the barrier modelEj 

There is, in addition to the non-uniform coupling to 
the magnetic field, an additional effect arising from the 
presence of crystallites with radius r < £(i,T). We shall 
show below that the magnitude of the jump in x"(w,t) 
will also depend upon ui (i.e. violate scaling) when the 
spin glass correlation length becomes of the order of, or 
larger than, the size of a spin glass particle. Conversely, 
the dependence of the change in upon u> over 

the full frequency regime could be used to generate the 
particle size distribution. 

The effective waiting time, tffi , after a magnetic field 
change, is given by Eq. (1), allowing one to write, 



,eff _ . ( E z \ 

t w -t w expy i 



(4) 



This scaling was first established by Vincent et aZ0 and 
by Chu et alB 

This scaling can be incorporated into the ac suscepti- 
bility. Before a dc magnetic field change, the relaxation 
of t) is well accounted for by a power lawjiil 



x>,t)-x:,ocM 



-b 



(•5) 



with b > 0, and x% the equilbrium susceptibility 
x"(w,*->oo). 

Using the relationship Eq. (4), the change in x"( u> it) 
upon a change in magnetic field at time t\ can be written 
as, 



Ax'WiHx'WfVx'Wi) 



(6) 



Here, tf^ is the effective waiting time upon a magnetic 
field change, defined through Eq. (4) with t w — > t%. 

With reference to Fig. 5, the lower the frequency 
oj, the greater the time t\ (because the abcissa is the 
scaling variable Lot). But the greater the time ti be- 
fore the magnetic field is changed, the larger the spin 
glass correlation length £ (t\ , T) , and therefore the more 
the likelihood that £(t%,T) > r, the size of the spin 
glass particle. But if £(ti,T) > r, the effective response 
time, tf, will be less than that for an "infinite" sample 
[£(ti,T) < r], and x"{ UJ i'^ ) w iH be larger, increasing 
the size of the jump from Eq. (6). However, ii) will 

also be slightly larger because it does not decay beyond 
x"(oJ,t erg ), where t erg is the time at which £(ti,T) = r, 
decreasing the size of the jump from Eq. (6). 

The increase in the first term in Eq. (6) turns out to be 
larger than the increase in the second because finite size is 



G 



involved in the argument of an exponential (through the 
Zeeman energy) in the first, while finite size enters only 
as an argument of a weak power law in the second (see the 
Appendix for details). This results in a larger magnitude 
of the jump in ^"(w,£i), the larger t\, or, equivalcntly, 
the smaller u>, precisely what is seen experimentally in 
Fig. 5. 

These arguments are qualitative, based upon Eqs. (4) 
and (5). A quantitative evaluation of Eq. (6) is made 
in the Appendix, fully supporting the conclusions of this 
subsection. 

Finally, for very short times, x"(o>,i) is seen to scale 
with ui-t in Fig. 5. Very short times means that £(t, T) < 
r for all particles. As we have already noted, under these 
conditions the barrier model calls for scaling with uj ■ t, 
again in agreement with experiment. 

V. KEY EXPERIMENTS AND NUMERICAL 
SIMULATIONS 

The arguments given above, especially in Sec. IV, are 
qualitative in nature, but fully capable of quantitative ap- 
plication. There are two approaches which we feel would 
be most relevant. The first, experimental, is one of care- 
ful magnetic field and waiting time variations upon a 
variety of samples. The second would be to make use of 
the model of Sec. II to numerically simulate particular 
realizations of spin glass materials. 

(a) Experimental determination of P(r) 

Preliminary examination of -a. number of spin glass 
samples, using SEM techniques,c3 suggest that many are 
made up of a powder-like array of small crystallites, em- 
bedding much larger apparently single crystal pieces. Of 
course, this division may not be general, and may only 
apply to those materials with which we have been work- 
ing. Nevertheless, the analysis of Sec. IV, parts (a) and 

(b) , can in this instance generate a measure of the bounds 
on the size distribution of the powder-like array compo- 
nent. 

The analysis of Sec. IV, part (a) generates the upper 
end, r max of the small crystallite length scale distribution 
P(r). The experiments are at fixed waiting time, with 
measurements as a function of the change in magnetic 
field, H. At very small magnetic field changes, only the 
volume occupied by £ 3 (t w ,T) contributes to E z , shifting 
the peak of S(t) by reducing the energy difference A — E z . 
As the magnetic field change increases, there will come a 
point when the energy difference A — E z becomes com- 
parable to A(r) — r max 3 XfcH 2 for the largest of the small 
crystallites. We have argued in Sec. IV, part (a), that 
this causes the "break" in the plot of logoff vs H 2 . 
Thus, the "break" field, Hb re ak, generates a determina- 
tion of r max . 

The other extreme of P(r), r m in, can be extracted from 
the departure from scaling, exhibited in Figs. 3 and 4. 



Here, the magnetic field change is fixed, and the experi- 
ments are a function of increasing waiting time, t w . From 
these two figures, the longer t w the more rapid the relax- 
ation of Mtrm (t, t w , T) as a function of the reduced time 

variable, — . For very small t w , £(t w ,T) is less than the 

tw 

"minimum" size of the powder-like array of small crys- 
tallites. The departure from scaling occurs first when t w 
increases to a point where £,(t w ,T) — r m i n , or equiva- 
lently A(t w ,T) = A(r min ). For longer t w , the crystallite 
with dimension r m i n is at equilibrium, and for that part 
of the sample, aging is over. Thus, the waiting time at 
which departure from scaling is first seen is a direct mea- 
sure of r m i n , the lower extreme of P(r). 

Careful (tedious!) measurements beginning from ei- 
ther domain, r max or r m ;„ , can of course be used to gen- 
erate all of P{r). At the very least, these two approaches 
will given a measure of the width of the P{r) distribu- 
tion. 

(b) Numerical Simulations 

An alternate, and certainly complementary approach, 
is to simulate the spin glass sample by selected choices 
for P(r). Previous simulations]!] using the barrier model 
of Sec. II, were able to duplicate the waiting time depen- 
dence of the response function S(t). Any attempt to fit 
to a particle size distribution P(r) will require the abil- 
ity to simulate logiotffl vs H 2 and MTRM(t 7 t w ,T) vs t. 
This procedure will not be unlike that of neutron or X-ray 
diffraction, where scattering from a specific model is mea- 
sured against experiment. It is the usual "inverse" prob- 
lem where small iterations from a hypothesized model are 
used to fit experiment. Previous success at fitting S(t) 
suggests that a similar procedure, using the data of Sec. 
Ill and the analysis of Sec. IV, will be successful. Having 
the limits r max and r m i n on P(r) in hand, as discussed in 
part (a) of this Section, will greatly aid such an analysis. 

VI. CONCLUSION 

This paper discusses spin glass dynamics for crystal- 
lites or amorphous particles of finite size. Departures 
from scaling arise when the spin glass correlation length 
becomes of the order of or larger than particle sizes. 
Qualitative arguments are given for the associated ex- 
istence of a "break field," flbreak away from a linear 

plot of logiQt^f vs H 2 ; the departure from — scaling 

tw 

of MTRM(t,t w ,T); and the frequency dependence of the 
magnitude of the jump in -x"{uj,t) vs u ■ t. A guide to 
future experiments and numerical simulations, leading 
to the extraction of the particle size distribution, P(r), 
are given with specific attention to what can be learned 
from experimental protocols. SEM measurementsES of 
the particle size distribution will lead to explicit experi- 
mental consequences, setting the stage for a consistency 
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check on the entire model. The authors find it remark- 
able that the behavior of magnetization measurements 
in the time domain could so directly depend upon the 
physical size parameters of the sample particulates. 

The authors have benefited from extensive discussions 
with Dr. J. -P. Bouchaud, experimental data supplied by 
Dr. M. Ocio, and from the financial support of the Japan 
Ministry of Education (Monbusho) and the U.S. National 
Science Foundation, Grant DMR 96 23195. 
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Appendix 

Sec. IV, subsection (c) gives a qualitative argument 
for the frequency dependence of the jump %"(w,i) vs tot 
upon a change in magnetic field as a consequence of finite 
spin glass particle size. This Appendix develops quanti- 
tative expressions for these quantities. 

To derive the dependence of the jump in x"(o>, t) upon 
change in magnetic field on particle size, consider two 
cases: I. r < £,(t,T) and II. r > £(i,T), where r is the 
radius of the spin glass particle, and £(£, T) the spin glass 
correlation length displayed in Eq. (2). We shall assume 
that the time dependence of x" is given by 

X " = X ': q (u) + A(u;t)- b , (Al) 

where x"ij( w ) m ay be different for cases I and II, but will 
cancel when we consider only the change in x" { UJ ^ t) upon 
a change in magnetic field. That is, we assume that the 
equilibrium (t — > oo) value of x"( w ) is magnetic field-hi- 
dependent. The exponent b « 0.20 from experiment.^ 

I. particle size r < £ (t, T) 

Consider the effect of a jump in magnetic field at a 
time t\, and for this subsection, assume that the particle 
size r < £(ti,T). Before the jump in magnetic field, 

XiMforei" ■ h) = X" QI H + A(u> • t erg r b , (A2) 
where we have used £(ii,T) = (ti/r ) a , with a = 
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0.153T/T g from Eq. (2), and where t erg = T r^ a is de- 
fined in Sec. IV, subsection (b). 

The effective time, tf after a jump in magnetic field, 
is given by 



,eff _. ( _ E z \ 

tj -t erg expy j 



(A3) 



This scaling was first-established by Vincent et a/0 and 
by Chu et at, givingEUI 



tfif 



T^r 1 1 a ex"p{^- 



knT 



(AA) 



Thus, after the jump in magnetic field, 

X'Uteri" ■ h) = X'UM + A{U ■ tfY h . (A5) 

Subtracting Eq. (A2) from Eq. (A5) gives the jump in 
x"(w,ti) upon a jump in magnetic field: 



Ax?(w,ti) = A(u>T r x ' a ) [earp( 
II. £(t,T) < particle size r 



br 3 Xfc H< 



(AS) 



Consider the effect of a jump in magnetic field at a time 
ti, and for this subsection, assume that £(ti, T) is smaller 
than the particle size. Before the jump in magnetic field, 



XuMfore^ ' *l) = Xe giI M + A(u • h)^ 

After the jump in magnetic field, 

Xll,after(u ' *l) = Xeq n M + • tjf)~ 6 

where, 
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fc R T 



(A7) 



(A8) 



(A9) 



The jump in X//( w ' *i) is then given by subtracting Eq. 
(A7) from Eq. (A8): 



the more the likelihood that £(ti,T) > r, the size of the 
spin glass particle. 

This means that the jump in x"{ u i *i) f° r case I should 
exceed the jump in x"(w,ti) for case II, or more simply, 
that Eq. (A6) should exceed Eq. (A10), more for smaller 
lu, or equivalently, larger t\. 

It is a somewhat tedious algebraic exercise, but one can 
show that this is indeed the case. Thus, finite size effects 
can generate the u dependence of the jump in x"i u> i 
This is a consequence of different E z values [through Eq. 
(A3)] as, a consequence of differing particle sizes. This 



featured adds to the non-uniform magnetic field cou- 
plings (larger, pthe larger the trap depth) introduced in 
the trap modelllil or (larger, the Wger the barrier height) 
introduced in the barrier model,E3 independent of possi- 
ble finite size effects, arising from the time dependence 
of T) and hence of E z . 



fba 3 r 3 Xf c H 2 
EXP { k B T 



(A1Q) 



For convenience, £ = ar, a > 1; a is a function of t\, T; 
and t\ — To(ar) 1 ' a . 

Experimentally, from Fig. 5, the magnitude of the 
jump in x"(w,<i) is larger, the smaller ui. But the lower 
the frequency u>, the greater the time t\ (because the 
abscissa is the scaling variable u>t). And the greater the 
time t\ before the magnetic field is changed, the larger 
the spin glass correlation length £(ti,T), and therefore 







